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Abstract. The magnetic properties of RM^Os multiferrroics as obtained by 
unpolarized and polarized neutron diffraction experiments are reviewed. We discuss 
the qualitative features of the magnetic phase diagram both in zero magnetic field 
and in field and analyze the commensurate magnetic structure and its coupling to an 
applied electric field. The origin of ferrolectricity is discussed based on calculations of 
the ferroelectric polarization predicted by different microscopic coupling mechanisms 
(exchange striction and cycloidal spin-orbit models). A minimal model containing a 
small set of parameters is also presented in order to understand the propagation of the 
magnetic structure along the c-direction. 
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1. Introduction 

The recent discovery of a new class of magnetic ferroelectric materials, in which electrical 
polarization coincides with a magnetic ordering or reordering transition, has generated 
very significant interest (El; The attractive feature of these materials is not 

so much the value of the electrical polarization, which is several orders of magnitude 
smaller than for typical ferroelectrics and even for "classic" magnetic multiferroics such 
as BiFeOa, but rather the very large cross-coupling between magnetic and electrical 
properties, which makes the "novel" multiferroics into enticing paradigms of functional 
behavior. Although conceptual devices based on novel multiferroics have been discussed 
(0), none of the materials so far described are directly suitable for applications, because 
the transition temperatures are low. What keeps much of the interest alive, however, 
is the possibility of discovering an underpinning general principle which could be 
applied to guide the synthesis of bulk of film materials with better properties. Some 
basic common facts about these materials have been established with clarity. As for 
conventional ferroelectrics, electrical polarization emerges as a result of a symmetry- 
breaking transition from a high-temperature paraelectric phase. Unlike conventional 
ferroelectrics however, the primary order parameter for this phase transition is magnetic 
rather than structural: as a result, the overall magneto-structural symmetry is lowered 
from that of the paramagnetic phase, eventually leading to a polar point group below 
one of the magnetic transition temperatures. Here, ferroelectricity is induced by some 
form of magneto-elastic and/or magneto-electronic interaction. On this principle, much 
work has been published on both group theory (J7|) and phenomenology (jg m, llO), 



aimed at establishing the symmetry requirements for the appearance of ferroelectricity, 
as well as the coupling between different order parameters and the generalized phase 
diagrams of these materials. Essentially all of these results are independent of the 
microscopic magneto-electric coupling mechanism. There is in fact no requirement 
that this mechanism is one and the same for all "novel" multiferroics — in this case, 
symmetry would be the single unifying principle. However, establishing this microscopic 
mechanism in each case is crucial for a quantitative understanding on the phenomenon, 
and it is an essential guiding principle for engineering new materials. As it happens, 
most of the "novel multiferroic" materials so far discovered share much more than 
a broken magnetic symmetry leading to a polar group: in fact, the presence of a 



cycloidal component to the magnetic structure in (lilt Il2t |4|; ll3|; ll5|) suggests a 



common underlying microscopic mechanism, which is critically dependent on the non 



collinearity of the spins. Nagaosa et al. ( 1l6l ) and M. Mostovoy ( llOl ) have explored this 



concept early on from different angles, and proposed what has now become known 
as the "theory of ferroelectricity in cycloidal magnets". In particular, Nagaosa et 
al. have proposed a detailed microscopic model based on the relativistic spin-orbit 
interaction, which is able to predict qualitatively, and to a certain extent quantitatively, 
the emergence of ferroelectricity in the presence of a cycloidal (or more generally non- 
collinear) magnetic structure of appropriate symmetry. It has been known for a rather 
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long time that other mechanisms can potentially induce ferroelectricity even in collinear 
antiferromagnets. However, the appeal of a single "universal" model is strong, and 
physics based on cycloidal modulations is very often sought as "the" single explanation 
for novel multiferroic behavior. The family of compounds with general formula RMn 2 05 
(R=Y, Rare earth, Bi and La) has so far stubbornly resisted being classified with the 
other compounds in the "cycloidal multiferroic" family. Three aspects of the RMn 2 05 
phenomenology stand out to suggest that these materials may be in a class of their 
own: the fact that the direction of the electrical polarization cannot be turned by an 



applied magnetic field (1171 ). the commensurate nature of the magnetic ferroelectric phase 
(ferroelectricity is all but lost in the low-temperature, incommensurate phase), and the 
fact that in the ferroelectric phase, moments in the afe-plane are almost collinear. The 
idea that RMn 2 Os may stand apart form the other multiferroics has been recently 
challenged by new findings — particularly the evidence for a previously unobserved 



cycloidal components (1181 ; Il9l ). In this paper, we describe the general phenomenology 
and present new results on the RMn 2 Os family multiferroics. We argue that, in spite of 
the emerging complexity and subtlety of the magnetic structure in these materials, there 
are still strong reasons to believe that cycloidal physics plays a minor role in RMn 2 Os 
(at least for the commensurate phase), and that ferroelectricity emerges due to more 
conventional exchange-striction effects in the context of a structure with built-in charge 
ordering. 

2. Experimental 

Poly crystalline, single-phase RMn 2 5 samples were prepared through conventional solid- 
state reaction in an oxygen environment. Stoichiometric quantities of Tb4C>7 purity 
99.998% , Dy 2 3 99.99% , Ho 2 3 99.995% , and Mn0 2 99.999% were thoroughly mixed, 
compressed into pellets, and then sintered at 1120 ° C for 40 h with intermediate 
grindings. The samples were finally cooled at 100 ° C per hour down to room 



temperature. Single crystals were prepared using the method described in (119l). Powder 
neutron diffraction data were collected on the GEM diffractometer at the ISIS facility 
of the Rutherford Appleton Laboratory (UK). The samples were enclosed in vanadium 
cans. All the data presented have been collected using either a helium cryostat or 
an Oxford Instrument 10 Tesla cryomagnet. For zero- field experiment, data have 
always been collected on warming, after cooling the sample to base temperature 
(typically 1.6K). Experiments performed in magnetic fields were conducted as follows: 
Meaurements on the YMn 2 Os compound were carried out at 1.6K in field of respectively 
0,2,4,6,8 Tesla. For each magnetic field, data were acquired for 2 hours. Measurements 
on the TbMn 2 C>5 system were carried out in the temperature range 20 to 34 K, in steps of 
IK, for several values of the magnetic field (0,1,3,5,7,9 Tesla). For each measurement, the 
sample was first zero-field cooled. The magnetic field was then applied and measurement 
were carried out on warming, setting a counting time of 1 hour per temperature. 
Single crystal neutron diffraction experiments have been recorded at the Institut Laue 
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Langevin (France), using the D10 4-circle diffractometer for unpolarized work and the 
D3 instrument equipped with CRYOPAD for Spherical Neutron Polarimetry. Details 
of the experimental procedures can be found in (Il9l ) and (120I ). Analysis of powder and 
single-crystal neutron diffraction data were carried out with the software FullProffj2ll). 



3. Crystal structure 



The crystal structure of RMn 2 Os compounds, already described in details in (1221 ; Il9l ). 
will be briefly reviewed here. The Mn ions are fully charge ordered, with Mn 3+ and 
Mn 4+ ions occupying sites of different symmetry. Mn 3+ and Mn 4+ are respectively 
coordinated by five oxygens in square pyramid geometry and six oxygens in octahedral 
geometry, as shown in Figure HJ The crystal structure is best described by considering 
the afe-plane configuration and out-of-plane configuration independently. In plane, 
octahedra and pyramids are corner-sharing through either the pyramid base or pyramid 
apex. In addition, adjacent pyramids are connected through their base. In total, three 
inequivalent exchange paths between magnetic ions exist in the plane (noted J3, J 4 and 
J 5 in Fig. [1] following notations of (0)). Along the c-axis, octahedral sites are sharing 
edges, forming linear chains. The Mn 4+ are located at z~0.25 and (l-z)~0.75, so that 
Mn 3+ ions, positioned at z=|, form layers in between adjacent Mn 4+ . R 3+ ions occupy 
sites of same symmetry than Mn 3+ in the z=0 plane, forming another layer alternating 
with the Mn 3+ layer. The ions in the primitive unit cell are labelled following the 
convention defined in (1221; Il9l). as indicated in Fig. [TJ 



4. Magnetic ordering -sequence of phase transitions 

The RMn^O ,^ st ructu re ty pe has been synthesized with Y, all the lanthanides (excluding 
Ce) (1231 ; [2J; |25| ). Bi (1261 ) and La (1231 ; 1271 ). Compounds containing lanthanum and all 
the rare earth lighter than Nd (included) do not become ferroelectric. The "typical" 
features of the magnetic phase diagram of the RMn 2 05 family are best described based 
on the temperature-dependent neutron diffraction patterns in Fig. [2j We can identify 
three or four regions depending whether the R site is magnetic or not. In the high- 
temperature paramagnetic phase, only nuclear Bragg peaks are present. On cooling, 
magnetic Bragg peaks first appear in reciprocal space positions corresponding to an 
incommensurate propagation vector q = (| — 8 X , 0, j + 5?). _ The yalues^ of 8 T . ( resp ^ S r ) 



29|; lafl; |3l|; l25|; |22|; |32j; |33j; M). 



range from 0.012 to 0.028 (resp. -0.013 to 0.027) (J2 
Within this high-temperature incommensurate phase(HT-ICP), the propagation vector 
changes with temperature, as it can clearly be seen from the curved shape of the 
Bragg peak centroid (this is best seen on the peak (0 1 0)±q around 7 A). On further 
cooling, the values of 5 X and 5 Z smoothly approach zero, and the magnetic structure 
"locks" into a commensurate phase (CP). The CP is the main ferroelectric phase for 
most compositions. Below ~ 20K, there is an abrupt transition to a low-temperature 
incommensurate structure (LT-ICP) — the position of the Bragg peaks continues to 
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change down to the lowest temperatures. Upon entering the LT-ICP, the value of the 
electrical polarization drops to a much smaller value (J3j), although there seems to be 
consensus that P ^ in some cases (e.g., Y). For Tb, as well as other magnetic rare 
earths, the intensity of the magnetic peaks rapidly increased below ~ 10K, indicating 
that the rare earth sublattice is becoming magnetically ordered (rare-earth ordered phase 
- ROP). The ROP is more strongly ferroelectric than the LT-ICP, and displays very 
large magneto-electric effects ( 35 ; [3); 0)- This general magnetic phase diagram, with 
some subtleties, is shared by RMn 2 5 with R=Y, Ho,Tb and Er (H; H; H; 0; S S), 
which are the most studied compounds. Excluding the obvious absence of rare earth 
ordering for Y and the somewhat different transition temperatures, the most significant 
differences between these compounds are within the HT-ICS. For Ho and Tb, the two 
components q x and q z of the propagation vector seem to become commensurate at the 
same temperature, whereas for Er and Y there is a small region where q z = \ while 
q x remains incommensurate (1371 ). The connection between this phenomenon and the 
appearance of "weak" ferroelectricity have not thoroughly been explored. In fact, hardly 
anything is know about the HT-ICS other than the propagation vectors, as the small 
magnetic moments make neutron studies more difficult than in the other phases. Three 
more compounds were studied in some detail: Non-ferroelectric LaMn 2 05 orders with a 
propagation vector q = (0, 0, ^) (1271 ; l33l ); DyMn 2 5 behaves like the other systems with 
magnetic rare earths above 8.4 K, but displays commensurate ordering with q = 0, 0) 



below this temperature, coexisting with a small fraction of the LT-ICP (j39t l33l ; 1341 ) ; the 



magnetic structure of BiMn 2 Os is commensurate and ferroelectric at all temperatures 
with a propagation vector q = (~, 0, ~) (1261 ; Il9l ). 



5. Influence of an applied magnetic field 

The strong magnetoelectric coupling in these materials can be directly evidenced by 
measurements of electric properties under magnetic field. For example in TbMn 2 05 
the upward jump in the dielectric constant at the CP to LT-ICP transition is pushed to 
higher temperatures as a magnetic field is applied (0), suggesting that the LT-ICP phase 
is stabilized by a magnetic field. This can be directly verified by constructing the H-T 
phase diagram of the CP-LT-ICP magnetic transition from powder neutron diffraction 
experiments in field. Data collected under fields up to 9 Tesla are reported in Fig. [3j 
Zero-field data show that the transition from CP to ICP phases is first order, as expected 
from group theory since the CP and ICP phases are characterized by wave-vectors 
belonging to different symmetry points of the Brillouin zone. The region of co-existence 
is very narrow and can be estimated to be about 1-2K. Under application of a magnetic 
field, the CP-ICP transition temperature is enhanced from ~ 25K in zero-field to more 
than 27K at 9 Tesla. The H-T phase diagram, reported in Fig. 4, has been constructed 
by fitting the data with two peaks of constrained widths, at positions fixed to the CM 
and ICP peak positions in zero-field. The transition temperature has been defined as 
the point of crossover between the intensity of these peaks. The stabilization of the LT- 
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ICP phase under magnetic field is in perfect agreement with the magnetic and dielectric 
phase diagrams (31), showing a gradual shift towards high temperatures, on one hand, of 
the kink observed in the first derivative of the magnetic susceptibility and, on the other 
hand, of the upward jump in the dielectric constant. The CP to LT-ICP transition is 
also smoothed out by magnetic field. Under field, the temperature range of the CP/LT- 
ICP phases co-existence is extended, while the propagation vectors of both phases do 
not seem to be altered. Overall, the magnetic behavior obtained by powder neutron 
diffraction data shows that the low temperature ICP phase is stabilized by application 
of an external field. This result is opposite to that recently published on the analog 
compound HoMn205 (l40l). showing that the boundary between CP and ICP phase is 
shifted to lower temperature under application of a magnetic field along the fe-axis. In 
both cases (Tb, Ho), neutron data correlates directly to the electric properties under 
magnetic field, confirming fundamentally different behaviors. Nevertheless, it should be 
noted that experiments reported for the Ho system (1401 ) and here have been conducted 
in different conditions. Work on H0M112O5 was performed on a single crystal and only a 
magnetic field parallel to the 6-axis was found to induce the ICP-CP transition. Work 
on powders is more qualitative due to the random orientation of crystallites with respect 
to the magnetic field direction. However, the general trend should not be affected by 
powder averaging, since magnetic fields along a and c proved to affect only slightly 
the magnetic state of TbM^Os . On the other hand, for HoM^Os, the induced 
CP state at high field is observed only in field cooled experiments, whereas zero-field 
cooled process leave the low-temperature ICP phase unchanged up to 13 Tesla. This 
indicates a large hysteresis (greater than 2 Tesla), as pointed out by Kimura et al., due 
to the first order nature of the transition. This does not explain, however, the opposite 
variation of the CP-ICP temperature boundary in the H-T phase diagram, and it is more 
likely that the discrepancy is a direct consequence of different single-ion anisotropies of 
the rare-earth ion (R) and its related effects on the Mn(d)- R(f) coupling. Another 
perspective in the comparison between Ho and Tb systems is given by their magnetic 
behavior at low temperatures. Tb ions order magnetically below 10K with the same 
propagation vector describing ordering of the Mn sublattice, i.e. k~(0. 48,0,0. 31) (|32l ). 
Under a moderately low magnetic field (2.5 Tesla), Tb orders ferromagnetically, i.e at 
k=0, without any noticeable change in the magnetic structure of the Mn ions. This 
behavior is in sharp contrast with the results on Ho (l40l ). where it appears that the 
magnetic ordering of the Mn and Ho sublattices arise at a single point in reciprocal 
space (CP or ICP) irrespective of the value of the magnetic field, suggesting a stronger 
d-f coupling. Moreover, the ferroelectric properties of TbM^Os, can be understood by 
the simple superimposition of several order parameters associated with the Mn and Tb 
sublattices respectively (Jl), implying the absence of coupling terms between them. We 
suggested that magnetic ordering of the Mn ions in TbM^Os is unaffected by the field- 
induced FM ordering of Tb based on indirect evidence, since the substraction of the 15 
K zero-field data from the 1.5 K, 2.5 T data, left an essentially perfect ferromagnetic 
pattern. This behavior can be verified by in-field neutron diffraction experiment on the 
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analog YMn 2 05. This compound shows the exact same sequence of CP-ICP transitions 
on cooling, but containing a non-magnetic R site. Figure [5] presents powder neutron 
diffraction patterns obtained in zero-field and in field up to 8 Tesla. The positions of 
the magnetic Bragg peaks remains unchanged under magnetic field, indicating that the 
propagation vector of the Mn magnetic structure does not vary. Also the magnetic 
peak intensities do not change within the error bars, confirming that the orientation 
of the magnetic moments remains unchanged. These results are in agreement with 
magnetization measurements in very high magnetic field, showing an extremely stable 
AFM structure at low temperature and absence of metamagnetic transitions for the 
Mn sublattice up to 40 Tesla. Taken together with previous experiments on TbMn 2 05, 
this tends to favor a picture in which the Mn and Tb magnetic sublattices are fully 
decoupled, since the metamagnetic transition of Tb does not influences the magnetic 
ordering of Mn. The situation for the Ho compound is therefore significantly different. 



6. Magnetic structures 



The magnetic structures of RMn 2 05 -type compounds have been studied since the 



seventies (14 It l42l ). and later in the eighties (1431 ) when the interest for them was mainly 



their unusual complexity. A continuous improvement of experimental and analytical 
techniques has lead to a refinement in their understanding, which, in the case of 
the incommensurate phases, has not reached completion yet (see section [TTT) . The 
magnetic structure is best understood in terms of magnetic a — b planes that are coupled 
along the c-axis. All the ferroelectric phases share closely-related in-plane structures, 
characterized by the doubling (for the CP) or near-doubling (for the ICP) of the 



magnetic a axis, while for non-ferroelectric LaMn 2 05 (propagation vector (0, 0, =) ( 1271 ) 
the a axis is not doubled. In particular, for all the commensurate, strongly ferroelectric 
phases, the in-plane magnetic structures are essentially identical. This suggest that the 
in-plane arrangement of the spins plays a key role in inducing ferroelectricity. On the 
contrary, the stacking of magnetic la yer s along the c axis varies considerably from one 



compound to the next — BiMn 2 Os (]26l : Il9l ) and DyMn 2 Os (1221 : l33t 1391 ) have a simple 



■ • ■ H 1 — ■ ■■ and • • ■ + + + + ■ ■■ staking, respectively, whereas for the other phases 

the stacking is either exactly (CP) or approximately (ICP) • • • + H h H 

As a typical example, the CP structure of YMn 2 05 is shown in Fig [6] — the figure 
also illustrated the symmetry-inequivalent magnetic interactions and their coupling 
constants. In understanding how this structure is established, the following factors 
need to be taken into account: 

(i) Exchange interactions. Although no detailed calculation has been presented so 
far, the in-plane interactions are likely to be uniformly antiferromagnetic. The 
interactions along the a direction are particularly strong, so that the presence of 
AFM zig-zag chains running along this direction is a common motif of all phases. 
For the LT-ICP, all magnetic structure determinations so far presented agree in 



8 



evidencing a long-wavelength modulation of this motif, although the details of this 
modulation differ. 

(ii) Frustration. The simultaneous presence of antiferromagnetic interactions and of 
fivefold exchange rings in the crystal structure immediately leads to a frustrated 
situation, where not all these interactions can be satisfied simultaneously. In the 
CP and ICP phases, this situation is relieved in different ways, possibly involving 
relaxation of the lattice and electronic structures and phasing of the zig-zag chains. 
Some of the resulting magnetic structure are acentric and capable of supporting 
ferroelectricity through magneto-elastic or magneto-electronic coupling. Although 
the matter is still controversial, we believe that this relaxation provides the primary 
explanation for the ferroelectric phase diagram of these materials. 

(iii) Anisotropy It is well-known that Mn 3+ ions with low-lying t^g orbitals tend to be 
magnetically anisotropic. For octahedral coordination (e.g., in LaMnOs, (44)), the 
anisotropy is related to the Jahn- Teller distortion — the spin direction parallel to 
the unfilled d z 2 orbitals is favored. In TbMnC>3 and other cycloidal multiferroics, 
the interplay between anisotropy (favoring collinear magnetic structures) and 
competing interactions (favoring cycloidal magnetic structures) is thought to be 
primarily responsible for the magnetic phase diagram — the former and latter 
arrangement are prevalent at high and low-temperature, respectively. In RMn 2 05 
, Mn 3+ has a pyramidal coordination, where the effect of the crystal field is very 
similar to the Jahn- Teller octahedral case. Consequently, the RMn 2 05 compounds 
also display significant anisotropy: in the CP, spin in each of the zig-zag chains are 
parallel to each other and to the axes of the pyramids. The axes of the pyramids are 
rotated by about ±15° away from the a axis, explaining why the a axis component 
of the spins, S x is much larger than the other two and also why the CP magnetic 
structure is slightly non-collinear. 

(iv) c-axis stacking. As shown in Fig. [1] and Fig. [BJ the stacking along the c axis is 
mediated by crystallographic chains of edge-sharing Mn 4+ 06 octahedra. To avoid 
confusion with the zig-zag chains running along the a axis, we will refer to these 
structures as "c axis ribbons". Here, there are two relevant magnetic interactions: 
through the Mn 3+ layers (J 2 ) and through the R layers (J\) as shown in Fig. [7J 
The Mn 3+ - Mn 4+ superexchange interactions through the common oxygen atoms 
appear to be the strongest contributors to J 2 which is therefore always ferromagnetic 
regardless of the sign of the in-plane interactions J 3 and J 4 (Fig. [[]). Consequently, 
in all phases, the spins in Mn 3+ ions adjacent through a Mn 4+ layer are always 
close to be parallel. We take advantage of this to simplify the description of 
the c— axis stacking, so that a single sign + or — actually represents two quasi- 
parallel Mn 3+ sites. The situation with the coupling J\ through the R layer is 
significantly more complex, explaining the diversity of the c-axis components q z of 
the magnetic propagation vectors. In the case of BiMn 2 Os and DyMn 2 Os , J\ is 
AFM or FM, respectively, and strong enough to enforce antiparallel ( 1 1 , 
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propagation vector q = (|,0, -)) or parallel (••• + + + + •••, propagation vector 
q = (|,0, 0)) staking, respectively. In all the other cases, J\ must be small, 
since both quasi-parallel and quasi-antiparallel stackings are found within the same 
magnetic structures. 

NNN and antisymmetric interactions At present, it is not known what determines 
the c— axis stacking, and the corresponding q z in materials other than DyMn 2 Os 
and BiMn205 . Two main effects are thought to be responsible for complex stacking 
and incommensurability in magnetic materials : competition between nearest- and 
next-nearest-neighbor interaction and the antisymmetric Dzyaloshinskii-Moriya 
(DM) coupling. The former is the accepted mechanism to explain the cycloidal 
phases in TbMnC>3 and other cycloidal multiferroics. An analogous model that can 
be applied to RMn 2 5 will be discussed in section |HJ It has recently been shown 



(I18l;ll9l) that in the CP of several RMn 2 Os compounds, including BiMn 2 05 , a small 
c-axis component S z of the spins is present, and that in the first approximation this 
is out-of-phase with respect to the other two components, giving rise to very flat 
cycloids running along the c axis. It is interesting to remark that this component 
is present regardless of the propagation vector, although it is smaller in the case of 
BiMn 2 05 . For BiMn 2 Os , there is a very simple explanation for the appearance 
of a c-axis component of this kind, due to the antisymmetric DM interaction (see 
Section [9] below) . It is tempting to extend this explanation to the other phases as 
well, and conclude that CP flat cycloids are of DM origin. We remark, however, 
that no detailed models or calculations have been presented so far. 

Leaving the small S z component aside, the general features of the in-plane magnetic 
structures of RMn 2 Os in the CP regime can be rationalized very well based on points 
(i) - (iii) above. From Fig. [1] and Fig. [6j it can clearly be seen that in the absence 
of structural distortions, there is an exact cancelation of all energy terms containing 
J3, since pairs of ions related by inversion symmetry carry magnetic configurations of 
opposite sign. This situation gives rise to the so-called magnetic Jahn- Teller effect, due 
to the analogy with the well-known structural effect. Here, the system can always gain 
energy from distortion of the crystal or electronic structure that makes the J3 slightly 
inequivalent, since the energy gain is linear and the energy cost is typically quadratic in 



the distortion. We have previously argued (1321 ; 1381 ) that this is the primary mechanism 



leading to the appearance for ferroelectricity in RMn 2 Os (see next section). 



7. Origin of ferroelectricity 

Two main models have been proposed to explain the appearance of ferroelectricity in 
the CP of RMn 2 05 ; here, we will refer to these models as the symmetric exchange- 
striction model and the cycloidal model, respectively. These two models are at present 
very difficult to disentangle based on the facts know to us with confidence. In particular, 
the point-group symmetry of the CP, which has been solved with great accuracy, is still 
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very high (m2m, see Ref (0)), so the ferroelectric polarization must lie along the b axis, 
regardless of the mechanism. 

The exchange-strict ion model relies on the magnetic Jahn- Teller effect as its main 
ingredient: here, electrical polarization would arise from a combination of atomic 
displacements and electronic rearrangements, removing the exact exchange degeneracy. 
The main features of the exchange-striction model are: 

• The most important element controlling ferroelectricity is the in-plane components 
of the magnetic structure. Non-collinearity is not an essential ingredient, as it has 
been shown that the same effect can be obtained in a collinear acentric structure. 
The stacking along the c axis is immaterial and so is the small c axis component, 
provided that all the layers have the same polarization. 

• Quantitatively, the electrical polarization is proportional to the scalar product of 
spins related by J3 in different zig-zag chains. Therefore, the relative phasing of 



the two chains is directly related to the magnitude of the polarization (138l ) . In 
the CP, this relative phase is uniform, so that each J3 pair contributes equally to 
ferroelectricity. However, there is no reason for this to be so in the LT-ICP — in 
fact for all models of the ICP so far proposed the scalar product varies in both 
sign and magnitude along the propagation direction. This provides an immediate 
explanation to the sudden loss of ferroelectricity in the LT-ICP (see below). 

As mentioned before, in the CP the electrical polarization is always directed along 
the b axis by symmetry. With all probability, there is a similar requirement for the 
incommensurate phases, so that whatever residual polarization remains in the LT- 
ICP should also be directed along the b axis, in agreement with all the experiments. 



Very recent electronic structure calculations have provided (1451 : 1461 ) evidence that 



ferroelecticity can be explained without invoking the spin-orbit coupling. 

The cycloidal model relies on the same spin-orbit-driven, inverse Dzyaloshinskii- 
Moriya effect that has been used to explain the physics of TbMnOs and other cycloidal 
magnets: here, electrical polarization would arise again from a combination of atomic 
displacements and electronic rearrangements, but this time what is minimized is the 
antisymmetric interaction between non-collinear spins. The general expression for the 
polarization thus generated is P oc e 12 x (Si x S 2 ), where Si and S 2 are spins on adjacent 
sites along the c axis and ei 2 is a unit vector connecting them (see for example (pi)). 

The main features of the cycloidal model are: 

• It provides a unified, "universal" explanation for all novel multiferroics. This is 
perhaps the most important reason of its popularity. 

• Non-collinearity in the context of a cycloidal structure is an essential ingredient of 
the model. 

• Locally, the polarization direction is defined by the cross product of two vectors, u, 
which is perpendicular to the plane of rotation of the spins, and q', defined as the 
projection of the propagation vector q along the direction of the bonds contributing 
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to antisymmetric exchange, q' may be difficult to define in general, because several 
bonds could provide independent contributions, but in the case of the CP, q' is 
directed along the c axis. The contribution to the electrical polarization of different 
ribbons may cancel in some directions. For example, in the CP, the plane of rotation 
of the flattened cycloids is defined by the c-axis and by the anisotropy direction, 
which is rotated 15° away from the a axis. However, the large P x components cancel 
out between the two ribbons in the unit cell, whereas the smaller P y component 
add, yielding an overall polarization along the b axis. Once again, this is entirely 
due to symmetry. 



• Preliminary investigations (1471) indicate that the LT-ICP has lower symmetry than 
the CP, but the polarization is still required to lie along the b direction. A LT-ICP 
precise solution of the magnetic structure may be able to discriminate conclusively 
between the two mechanisms. 

Most of the consideration made so far are qualitative. However, it is possible to 
make a quantitative assessment of the likely contributions of both exchange-striction 
and cycloidal mechanisms, once the magnetic structures are known in detail. What 
can never be known from the magnetic structure alone are the coupling constants for 
symmetric and antisymmetric exchange, although it is reasonable to assume that the 
antisymmetric coupling should be weaker, since it relies on the intrinsically weak spin- 
orbit mechanism. In both cases, the electrical polarization P is obtained by multiplying 
the appropriate coupling constant times a polar vector (E or S below), which is entirely 
defined by the magnetic structure and has the dimension of fi 2 B . For the calculation of 
these polar vectors, only the spins within one unit cell need to be taken into account, 
once the propagation vector is known. Therefore, the same formalism can be applied to 
the CP and to all the ICPs. For the calculations presented in this section, we employ 
the following convention: the three components of the magnetic moment of a given site 
s, Vf (i=x,y,z) in unit cell R; are: 

V:(R l )=M l s cos(qR l + ( p s i ) (1) 



The MjS and <piS can be calculated from the Fourier coefficient given in (1191 ) . We note 
that the pro pag ation vector q and phase <fi are given here in radians and not in fractional 
units of 2ir (19). The same labels than H 0) are employed for magnetic ions in the 
unit-cell. 

The exchange striction polar vector E y is always directed along the b axis. If q x is 
incommensurate, E y is calculated as 



/ i=l,3 

- Mf Mf cos(0f - 0f + q x ) - Mf Mf cos(# - <f>f + q x ) 
+ Mf Mf cos(0f - <j>f) + Mf Mf 008(0}* - 0f ) 

- Mf Mf cos(0f - 0f) - Mf Mf cos(0f - <f>f ) (2) 
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If q x = ir, there is an additional "umklapp" term, so that 



E y = J2 M f M t cos <pf cos <j>f + Mf Mf 1 cos <j>f cos . 

i=l,3 



- Mf Mf 4 cos0f cos(0f + & ) - Mf Mf 4 cos0f cos(0f + q x ) 
+ Mf Mf cos0f cos 0f + Mf Mf cos 0f cos <f>f 

- Mf Mf cos 4>f cos 0f - Mf Mf cos 0f cos <j>f (3) 

There are two separate components of the spin-orbit polar vector, referring to the 
cross product of spins across the Mn layer (S 1 ) or across the RE layer (S 2 ). The same 
formula serves to calculate both x and y components. For q z 7^ 7r we obtain: 



1 1 

1 2 


Mf Mf 


cos( 




0f) 


- Mf Mf cos( 


0f- 


€) 


+Mf Mf 


cos(( 


r - < 




- Mf Mf cos(^ 


6f _ 


€) 


2 1 

* 2 


Mf Mf 


cos( 


*?- 




- fo ) - Mf Mf 


cos(( 


pf- 


+Mf Mf 


cos(( 


t - < 


/f - 


- 9z )-Mf Mf 


cos((/ 


f - 
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Whereas for the umklapp case q z = ir of BiMn 2 05 we obtain: 



(4) 



5/ = Mf Mf cos 4>f cos 0f - Mf Mf cos 0f cos </>f 

+ Mf Mf cos 0f cos <f>? - Mf Mf cos <j>* cos 0f 
S? = S} (5) 

The total polarization is given by 

P _ r .SO ql _|_ SO o2 

p 2/ = c es ^+ c r^+^°^ 2 (6) 

where cf , and c es are magneto-elastic coupling constants and can have either 

sign. 

Calculated data for the different polar vector components based on published 
magnetic structures are reported in Tab. [TJ The following observations can be made: 

(i) In the strongly ferroelectric CPs, the exchange-striction polar vector E y is always 
much larger than the cycloidal polar vectors — typically by a factor of 50-100. 
This is in itself a strong indication that exchange striction should be the dominant 
mechanism if one considers the fact that the exchange-striction coupling constants 
are also expected to be larger. 
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Compound q E y S% Si S% Sy 



YMn 2 5 CP (19) 


(|,o, J) 


22.6 


-0.01 


-0.16 


-0.04 


-0.39 


YMn 2 5 ICP (48) 


(0.48,0,0.29) 


6.8 


-5.84 


1.34 


-4.88 


-0.49 


HoMn 2 5 CP (19) 




27.4 


-0.10 


-0.23 


-0.21 


-0.46 


BiMn 2 5 CP (19) 


(I,o,|) 


-45.4 


-1.38 


0.08 


-1.38 


0.08 



Table 1. Magnetic polar vector components relevant for the exchange-striction 
mechanism (E y ) and the spin-orbit mechanism S y , S% and Sy), calculated from 
Eq. OSfor different RM^Os compounds. We employed our data for the CPs, whereas 
for the ICP of YMn 2 05 we used recent data by Kim et al. ifiih. All values are in /j? B . 
The different signs of E y are not significant, since they are due to the inversion-domain 
variant chosen for the refinement. 



(ii) We have reported for comparison the polar vector components for the most recently 



published refinement of the YMn 2 Os ICP by Kim et al. (1481) . This structure 
contains cycloids in the a — b plane as well, which could in principle contribute 
to the ferroelectricity through the spin-orbit coupling mechanism. We observe, 
however, that the strong decrease of E y is in accord with the drop in the electrical 
polarization, as observed experimentally. The decrease of E y is due to the fact that 
the zig-zag chains are now out-of-phase rather than in-phase, as we remarked for 



the ICP structure we previously obtained from powder data (138l ) . On the contrary, 
the trend on the cycloidal polar vector components — in particular the prediction 
of a significant polarization along the a axis — is inconsistent with the experiments. 

(iii) The large value of E y calculated for BiMn 2 Os is in agreement with the 
experimentally observed large polarization. However, our magnetic structure would 
be consistent with a spin-orbit polarization directed along the x axis, at odds with 
the experiments. Nevertheless, a more precise determination of the s z components, 
which are very small for BiMn 2 05 , would be required to be able to employ this 
argument with confidence to validate the exchange-striction model. 

8. A simple model for the propagation along the c axis 

The magnetic phase diagram of the RMn 2 05 compounds has been studied 



phenomenologically in great detail by Harris and coworkers ( 149l ; l50l ). In this section, we 
propose a more modest approach, based on the well-known linear-chain model, which, 
however, contain the minimal ingredients — competition between nearest-neighbor (NN) 
and next-to-nearest neighbor (NNN) interactions, widely thought to be responsible for 
the observed cycloidal structures in TbMn03 and related compounds. The purpose of 
this discussion is to ascertain whether these ingredients are sufficient to explain the 
general tendency of RMn 2 Os to develop an out-of-phase s z component for a variety of 
propagation vectors. 

The spins and interactions to be considered , as shown in Fig [Tl are: 
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s' 3 = ls x cos {j3R z - a/2 + (3) + ks z sin ((3 R z - a/2 + (3) 
s a = i s x cos (/3 i? 2 + a/2) + ks z sin (/3 _R 2 + a/2) 
82 = 15;,. cos {(3 R z — a/2) + k s z sin {(3 R z — a/2) 

s 3 = isa, cos (13 R z + a/2 - /?) + ks 2 sin(/?# 2 + a/2 - (3) (7) 

where i and k are unit vector in the xy-plane and z directions, respectively, R z 
is the unit cell index, j3 is the propagation vector component along z and a/2 is the 
phase of the spins, the latter two being expressed in radians, so that they relate to the 
tabulated values as a = 4ti5 and (3 = 2nq z (jig ). Here, s x and s z refer to a generic 
in-plane component and to an out-of-plane component, respectively. 

We can construct an idealized 2-dimensional model, coupled through NN 
interactions (J\ and J 2 , which are different) and NNN interactions (J'). The total 
energy per unit cell is 



E = J 2 Si ■ S 2 + ^ J\ (Si ■ S3 + s 2 • S3) + ^ J' (si • S3 + s 2 • s 3 ) 

Rz " ~ 

+ Ts 2 z (sin 2 ((3 R z + a/2) + sin 2 (J3 R z - a/2)) (8) 

The last term containing T represents the contribution of the anisotropy. By 
replacing the expressions of the spins from eq. [7] and decomposing the energy into 
the normal and umklapp terms E = En + Ey we obtain: 



s 2 + s 2 

E N = x - z ( J 2 cos a + Ji cos(a - (3) + J' cos (3) + Ts 2 z (9) 
2 2 

£c/ = E cos(2/5i? 2 ) ( J 2 + J x cos/3 + J' cos [a -13)) 

Rz 

s 2 

-r-^^cos(2/3 J R 2 )cosa (10) 

N Rz 

One can easily see from eq. [9] that the umklapp term is non-zero only for (3 = 
and (3 = n (doubling of the unit cell along the z axis). 

Let us first consider the general case in which the umklapp term is zero. We note 
that we can ignore the anisotropy when we minimize the energy as a function of the 
two angles a and (3, since the anisotropy term does not contain the angles. By defining 
R = J' I J 2 and T = J\j J 2 and setting the first derivatives to zero we obtain the two 
solutions: 



sin a 
with 
sin (3 



= — R sin (3 
= 
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or 

The first solution always gives rise to umklapp since (3 = or ir, so it should 
always be discarded. We can always set (3 > 0, whereas a is set in the appropriate 
quadrant to satisfy eq. [TTJ Once this is done, one should obtain the values of s x and 
s z under the constraint that the total spin at any site cannot exceed the full ionic value 
(3/2/xb for Mn 4+ ). For an elliptical cycloid of this type, the maximum spin value is 
max(s x ,s x ) = s x , i.e., the long semi-axis of the ellipse. Therefore, regardless of the 
specific values of a and (3, the minimum energy is always attained for either s z = s x 
(small anisotropy, circular cycloid) or s z = (large anisotropy, collinear spin density 
wave). 

The two situations riving rise to umklapp are (3 = and (3 = ir. We note that 
in this simple model we cannot obtain the (3 = n/2 lock-in situation that characterizes 
the ferroelectric phase. In order for this to occur, one would need a quartic term in the 
free energy — for example, through magneto-elastic interaction. The presence of these 
quartic terms manifests through the appearance of charge peaks at twice the magnetic 



propagation vector (15 If ). Setting p = s z /s x , the two expressions to minimize in the 



presence of umklapp become: 



E 



S: (Ji + J 2 )(l-P 2 ) + J'(l+p 2 ) + 2Tp 2 



2 

+ cos a 



(Ji + J 2 )(l + P 2 ) + J'{l-p 2 )-2Tp 2 



for (3 = (12) 



and 



E 



-Ji + </ 2 )(l-p 2 ) 
-Jx + J 2 )(l 



J'(l + p 2 ) + 2Tp 2 



P 2 W(i-p 2 )-2r P 2 



2 

+ cos a 

for (3 = 7T (13) 

With the help of eq. [9lfm [T2l and [13] one can construct phase diagrams for what is 
effectively a two-parameter problem. The case relevant for RMn205 is that of J 2 < 0, 
since, as we mentioned previously, the interaction through the Mn 3+ layer is always 
ferromagnetic. Fig. [9] shows the phase diagram at fixed J 2 = — 1 as a function of J\ 
and J' and in the isotropic case (s x = ±s z , T = 0). We have plotted only the region in 
which J 2 is the largest exchange constant — a plausible assumption, as we have seen, 
but the extended phase diagram is equally easy to plot. 

The following observations can be made by inspecting the phase diagrams: 

(i) Even in the absence of umklapp, there are two large regions of commensurability. 
For (3 = the magnetic cell coincides with the chemical cell, whereas for (3 = it the 
magnetic cell is doubled. 
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(ii) More significantly, a = always whenever (3 = or (3 = it. This is also clear by 
inspecting eq. [12] and [131, since this expression is minimized for one of the extrema 
of cos a. Since a = in this part of the phase diagram, the solution is a collinear 
AFM, and the terms containing p and T in eq. [T2l and [TBI cancel out. This is very 
important, since it indicates that this simple model is incapable of reproducing the 
BiMn 2 05 situation, where the presence of a s z component indicates that a > 
even with (3 — ir (propagation vector q = (|, 0, |). 

(iii) As soon as one departs from the isotropic situation, the umklapp terms tend to 
stabilize the commensurate phases, and the region of incommensurability shrinks. 
For p = 0.5, the portion of the phase diagram shown in Fig. [9] is completely 
commensurate with a = 0, and is split diagonally (top left to bottom right) between 
(3 = (bottom) and (3 = ir (top). Typical experimental values of p are even smaller 
— of the order of 0.2. 

(iv) Even in the presence of umklapp, one can always find small patches of 
incommensurability in the extended phase diagram. However, in order to 
stabilize the RMn 2 Os incommensurate phase by this mechanism, one would require 
extremely fine-tuning of the exchange parameters. 

One also observes that there is complete degeneracy between the configurations 
differing by the sign of the z components of the spins (a — > —a , (3 — > — f3), corresponding 
to the counter-clockwise and clockwise rotation of the cycloids (if one imagines to 
travel along the positive z direction, and to the two different polarities of the magnetic 
structure (there is no true chirality here). This is completely trivial in this simple ID 
model, but it is not so in the realistic situation of RMn 2 05 , since, as we have seen, 
the in-plane magnetic structure also has a polarity, and reversing all the z components 
would correspond to a non-symmetry equivalent domain (see Section [10] below). The 
degeneracy between these domains must be lifted to some degree, in order to enable the 
observation of a s z component by diffraction. 

This observation, together with those at point [u] and [Tv] above, suggest that an 
additional mechanism, in addition to competition between NN and NNN interactions, 
may be required to explain the observed phase diagram of RMn 2 Os . 

9. Dzyaloshinskii-Moriya interaction and its consequences for the RMn 2 Os 
magnetic structures 

As we have seen, the simple model presented in Section [8] is incapable of reproducing 
some important features found in RMn 2 Os compounds — for example the presence 
of a small s z component in BiMn 2 Os . Moreover, opposite cycloidal polarities (i.e., 
opposite signs of s z , everything else being equal) are energetically degenerate, while this 
degeneracy must be lifted to be consistent with the diffraction observations. In fact, this 
degeneracy is characteristic of all models that contain symmetric exchange terms only, 
at least for phases that retain part of the crystallographic symmetry. In fact, symmetric 
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coupling energy terms that split the degeneracy must be of the form S % x or S l y S J z , 
where i and j refer to different magnetic sites. If a set of mirror planes perpendicular 
to the c axis is retained, as it is the case, for example, for the commensurate phase of 
RMn 2 05 (1191 ). terms of this form cancel out exactly by symmetry. It is therefore useful 
to look specifically for an antisymmetric interaction term that is capable of lifting the 
degeneracy. 

In the crystal structure of the paramagnetic phases, the DM vector D between 
Mn 4+ atoms chained along the z axis is identically zero, because these sites are related 
by inversion (Fig. [8]). In the ferroelectric phase, this is no longer strictly true, because 
the ferroelectric displacements break inversion symmetry. However, these displacements 
are extremely small, and so should also be the resulting D vector. Here we show that 
there is a much more efficient mechanism to generate non-collinearity along the z axis, 
through an effective antisymmetric interaction mediated by the Mn 3+ atoms. Mn 3+ and 
Mn 4+ atoms are not related by any symmetry, and the D vector associated with pairs of 
such atoms is in no way restricted. However, D transforms like an axial vector between 
symmetry-related bonds, so some of the DM interaction terms cancel out by symmetry. 
The labeling scheme we employ to calculate the DM energy is shown in Fig. [8], and is 
consistent with (1191 ). 
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where D" and D are the components of the DM vector parallel/perpendicular to 
the ab plane. By employing Eq. [THwe can calculate the DM energy E DM : 



Edm — Dx_ 4 • [(So4 — S a i) x (S&i — S b2 )] + 
D^_ 4 • [(S a4 + S 0l ) x (S w + S 62 )] + 
D!_ 2 • [(S a2 - S a3 ) x (S 61 - S 62 )] + 
D^_ 2 • [(S a2 + S a3 ) x (S M + S w )] = 
~ 2 d|_ 4 ■ [S a4 x (S w - S 62 )] + 2 D^_ 2 ■ [S a2 x (S 61 + S 62 )] (15) 

In the last line of Eq. [TBI we have exploited the fact that S a4 ~ — S a i and S a2 ~ S a 3 
(the equality is exact for fully constrained models and q x = 1/2). One should also 
remark that for the commensurate phases S a2 is essentially in-plane — it is in fact 
exactly so by symmetry for BiMn 2 0s and approximately for most of the other elements 
with q = (i, 0, |). Therefore, the effect of the first term is clearly that of inducing a z 
component of opposite sign, and therefore a canting angle between S^ and S& 2 - It is 
an effective antisymmetric interaction along the chain, mediated by the in-plane spins. 
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The second term induces an in-plane canting angle S a and S&, spins. Since this does 
not break the symmetry, we will not consider it in the remainder. 

We can reproduce the effect of the first term in our simple 1-dimensional model 
by making a correspondence between spins in Fig. [7| and Fig. [8] and setting the DM 
vector along the b axis, i.e., d|_ 4 = }D. We also assume S a 4 = s x cos (3R Z . After some 
manipulations we obtain: 

2D 

Edm = -jT- s * s z cos P R z [sin {(3R Z - a/2) - sin ((3R Z + a/2)] = 
Rz 

= — Ds x s z sin a/2 (16) 

It is clear that Edm stabilizes a non-zero value of a even for the umklapp cases (Eq. 
IT21 and [T3l) . This simple mechanism is therefore capable of explaining the observation of 
a small s z component in the case of BiM^Os , where (3 = tt. We can see now that Edm 
discriminates between right-handed and left-handed cycloids, since its sign is reversed 
if one changes the sign of a. In addition, the sign of a (reversal of cycloidal polarity) 
must change in response to a change in the sign of s x . For a given s x , changing the 
sign of s x corresponds to reversing the in-plane polarity of the spin system. Therefore, 
the DM term favors a consistent alignment of in-plane and cycloidal polarities, always 
parallel or always antiparallel depending on the sign of D. This lifts the degeneracy 
between the two types of non-symmetry-equivalent domains — an essential ingredient 
to explain why cycloids are at all observed by diffraction. 

The effect of Eidm on the wider phase diagram, away from the umklapp points, 
is more difficult to calculate analytically. Preliminary numerical calculations suggest 
that, in addition to making a non-zero for the umklapp phases, the DM term has 
the additional effect of extending the region of incommensurability — unsurprisingly 
perhaps, because it promotes non-collinearity. 

10. Magnetic Domain switching under an applied electric field, Spherical 
Neutron Polarimetry 

The simple model outlined in Sections [8] and [9] paints a picture of the RMn 2 05 physics 
that is rather different from that of the typical cycloidal multiferroics. The case of 
BiMn2C>5 demonstrates with particular clarity that the cycloidal component does not 
emerge independently as for TbMnC>3 , but is induced by and has consistent polarity 
with the in-plane magnetic structure. The presence of an antisymmetric energy term 
that lifts the degeneracy between configurations with parallel and antiparallel polarities 
is also important, because it has a direct implication on the domain formation and on 
the domain switching upon application of an electric field. In the presence of this term, 
if the configuration (a, j3) is the stable state, the "alternate" configuration with reversed 
cycloids, (—a, —f3) is not an extremum of the energy function. However, if Edm is small, 
in general there will be a local minimum near the (—a, —0) point, which corresponds 
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to a metastable state. In this case, it is possible in principle to reverse one of the two 
polarities (in-plane or cycloidal) without affecting the other. Which of the two polarities 
is reversed depends on the specific coupling mechanism to the electric field; in other 
words, if the cycloids were entirely responsible for the electrical polarization through 
spin-orbit coupling, the direction of rotation of the cycloids must be reversed when the 
electrical polarization is reversed. Likewise, reversal of the in-plane polarity is a strict 
requirement of the exchange-striction model. For certain values of the parameter, the 
metastable minimum is very shallow or absent: in this case, we would expect that the 
electric field would switch the structure between truly degenerate domains related by 
inversion. We can summarize the different scenarios as follows (metastable states can 
be "shallow" or "deep" with respect to ksT): 

(i) FE by spin-orbit coupling - deep metastable state. The prediction for this scenario 
is that only the cycloids will switch upon reversal of the electric polarization by an 
applied external field. The in-plane structure should remain unaffected. 

(ii) Both spin-orbit and exchange-striction contribute to FE, deep metastable state. Here, 
both in-plane and cycloidal components should switch, possibly at different fields. 
This scenario is unlikely, as it should produce a 2-step FE hysteresis loop that is 
not observed experimentally. 

(iii) FE by exchange-striction, deep metastable state. Here, only the in-plane structure 
should switch, whereas the cycloids should not be affected. 

(iv) Any mechanism, shallow metastable states. Here, the two components always 
switch simultaneously, so only domains related by inversion symmetry are ever 
observed. 

A direct observation of the response of the magnetic structure to the reversal of 
an external electric field is therefore non-trivial and potentially informative about the 
mechanism of ferroelectricity. 

Very few techniques are available to probe directly the domain structure of an 
antiferromagnet. Second-harmonic light generation has been used successfully to this 



effect in multiferroics — see for example (1521 ) . However, this technique can only be 



used for "r-point" antiferromagnets, in which the magnetic unit cell coincides with the 
primitive chemical cell. In fact, the macroscopic quantity that is relevant for magnetic 
SHG is a third-rank, time-reversal-odd axial tensor. For non r-point antiferromagnets, 
time-reversed translations are symmetry operators. Therefore, the magnetic point group 
contains time reversal, and the tensor is identically zero. Another way of saying the 
same thing is to observe that SHG is an optical technique, and should therefore only be 
sensitive to zone-center (r-point) effects. In these more complex cases, the technique of 
choice to probe antiferromagnetic domains is scattering of polarized neutrons, of which 
Spherical Neutron Polarimetry (SNP) is a particularly powerful version. The neutron 
spin polarization is defined as the statistical average over the neutron beam of the 
expectation value of the quantum mechanical spin projection operator. An almost fully 
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polarised (> 99 %) monochromatic neutron beam is generated by means of a Heusler 
crystal monochromator. The neutron spin direction can be altered in one of three ways: 

(i) adiabatic rotation : the spins will follow a slowly-varying magnetic field; 

(ii) precession: when crossing an abrupt step change in the field direction, the neutron 
spins will precess around the new magnetic field; 

(iii) flip: the neutron spins can be flipped by 180° by a special device, known as a flipper, 
which also exploits neutron spin precession. 

By an appropriate combination of adiabatic rotations, precessions and flips, one can 
prepare the incident neutron polarization in an arbitrary direction, and also rotate the 
scattered polarization so that any of its components is parallel (or antiparallel) to the 
analysis direction. To this end, a flipper and a polarization analyzer, in this case a 3 He 
spin-filter, are placed in the scattered beam in front of the neutron detector. The sample 
itself is held in exactly zero magnetic field. Unlike conventional neutron crystallography, 
spherical polarimetry does not rely on measured neutron intensities, but on intensity 
ratios between parallel and antiparallel settings of the final flipper /analyser pair. In 
practice, one can measure 9 independent "flipping ratios", by setting the incident and 
scattered polarization along the X, Y and Z directions (in appropriate coordinates, see 



53)). All the elements of the resulting polarization matrix (for complete treatment see 



(1531 ; \54i)) are sensitive to the magnetic structure factor, as is the unpolarized neutron 



scattering cross section. However, some of the matrix elements - - in our case the 
elements P zx and P yx — carry unique information about the domain population, even 
if these domains are related by inversion. 

We performed a SPN experiment on a YM^Os single crystal as a function of 
temperature and external electric field applied along the h crystallographic direction. 
The results of this experiment are described in details in Ref. (120l). With the NPS 
technique, we can probe directly both domain structure and domain population for 
different temperatures and orientations of the external electric field. However, on D3, 
there is only limited access to reflections out of the horizontal scattering plane, so the 
crystal mounting will dictate which region of reciprocal space is observed. For this first 
experiment, we chose to mount the crystal with the b axis vertical, i.e., perpendicular 
to the scattering plane, so that only reflections with k = were accessible. This setting 
enabled to probe in great detail the magnetic structure within the a — b plane, and to 
observe the effect of the domain switching on the in-plane structure. We were, however, 
not able to distinguish between stable and "alternate" domains. Consequently, our 
experiment can completely corroborate or falsify scenario (i) above. It can also provide 
partial information about (ii), but it is unable to distinguish conclusively between (iii) 
and (iv). The main results of our experiment are summarized in Fig. [101 We applied 
two types of electric field-switching protocols: cooling in an applied electric field of 
either polarity (±2.2kV/cm ) through the Neel temperature T^r down to 25 K (left 
panel) and cycling the electric field at 35 K, i.e., closer to T/v, after polarizing the 
sample at 25 K. For in-field cooling, we observed a complete reversal of the in-plane 
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domain population (Fig. [TUl left). This clearly demonstrates that the in-plane magnetic 
structure is coupled to the electric field, ruling out scenario (i) above. At 35 K, we were 
able to measure a complete hysteresis loop on the off-diagonal neutron polarization 
matrix elements (Fig. [TUJ right). The loop is shifted downwards with respect to the 
centerline, indicating that some domains are "locked" in a fixed polarization direction 
and cannot be reversed by the small electric field available to us. More importantly, the 
SNP results are in excellent agreement with measurements of the macroscopic electrical 
polarization performed with the same protocol, including the downward shift of the 
hysteresis loop and limiting values of the electrical polarization at the top and bottom 
of the hysteresis loop. In other words, the in-plane antiferromagnetic domain population 
is strictly proportional to the electrical polarization. This result indicates that, if both 
in-plain and cycloids contributed independently to the polarization (scenario (ii)), the 
latter contribution is probably very small. In a future experiment we are planning to 
perform an independent measurement of the cycloidal domain switching under similar 
conditions, by mounting the crystal with a different orientation. This will enable us to 
discriminate among the above scenarios in a unique way. 



11. Outlook and conclusions 



As we have seen in the previous sections, much progress has been made towards 
understanding the connections between magnetism and ferroelectricity in RMj^Os , 
thanks to the sustained effort of several groups. Neutron diffraction played a major 
part in this research, providing several key pieces of information. Nonetheless, much 
remains to be done. The most obvious gap is the lack of a definitive determination 
of the magnetic structures of the two ICPs. The HT-ICP is clearly the most difficult 
one, because the magnetic moments are small, and no solution has been proposed to 



date. Starting form the seventies, several grou 



structure of the LT-ICP, from both powder (1421 : l38l ) and single-crystal data ( 1391 : l43l : 1481 ). 



3S have attempted to solve the magnetic 



Some features are common to all these refinements, and can be considered as established 
with good confidence. Other aspects are slowly emerging as techniques and refinements 
improve. The main results to date can be summarized as follows: 

(i) In the LT-ICP, the modulations of the two zig-zag chains are approximately in 
quadrature, rather than in phase, as in the CP. This emerges clearly both from the 



powder data (38) and the more recent single- crystal data (|48j). This is a key feature 
in the context of the exchange-striction model, since it explains why ferroelectricity 
is suppressed in the LT-ICP: in fact, all terms in exchange-striction polar vector 
E y (eq. EJ) contain the cosines of the phase differences between atoms on different 
zig-zag chains. If the phase difference between these chains is close to 90°, the 
cosines will be small and so will be the electrical polarization. 

In the LT-ICP magnetic structure, one can recognize "patches" resembling CP 



domains (38). Both polarities are equally represented, so their contribution to 
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ferroelectricity cancel out. These "patches" alternate with regions where the dot 
product between sites on the two chains is very small or zero, yielding a negligible 
contribution to ferroelectricity. 



(iii) The main difference between the powder structure 0381 ) and the more recent (and 
accurate) single-crystal structure |48l ) is the relative phase of the b component of 
the spins. In the latter, the S x and S y components are in quadrature yielding a 
cycloidal component in the a — b plane. Here, the spins of the two chains are 
almost orthogonal throughout the modulated strucure. This has no impact on the 
exchange-striction polar vector calculation, since the two components contribute 
independently to E y , but may provide an additional spin-orbit contribution. We 
remark that if the two zig-zag chain modulation were identical except for a phase 
factor — a reasonable hypothesis — this new spin-orbit contribution would be 
directed entirely along the b axis. 

(iv) Assessing the the spin-orbit polar vector in the LT-ICP requires a level of accuracy 
in the magnetic structure determination that, in our opinion, has not yet been 
attained. In Section [7J we present the results based on the recent determination 



by Kim et al. (]48l) . which are not in accord with the experimental values of the 
electrical polarization. We remark, however, that the structure by Kim et al. is not 
yet completely satisfactory, since symmetry-equivalent atoms have widely different 
modulation amplitudes. We believe that these inconsistencies should be resolved 
before a final assessment of the relative importance of the two mechanisms can be 
made. 

Another aspect so far remains largely unexplored is that of rare earth magnetic 
ordering. It is now known that a sizeable magnetic moment is induced on the RE 



by the ordering of Mn (see, for example, (I5ll)). Some aspect of the low-temperature 



ordering and the effect of the application of an external magnetic field have also been 
established ([32). However, the dramatic changes of the electrical polarization upon RE 
magnetic ordering have not been explained. The study of RE with different on-site 
anisotropy may provide important clues to understand these effects. 

In summary, we have presented an overview of the neutron diffraction results 
on the RMi^Os multiferroics, with particular emphasis on the correlations between 
magnetic ordering and ferroelectricity. Solely on the basis of the refined magnetic 
structures, we calculate polar-vector quantities that can be directly related to the 
electrical polarization, for both exchange-striction and spin-orbit microscopic models 
of the dominant magneto-electric interactions. We also explored the underlying causes 
giving rise to the complex magnetic phase diagram of RM^Os , and proposed a 
"minimal" model that can reproduce some of the observed features, whilst highlighting 
the importance of antisymmetric exchange in stabilizing phases of consistent polarity. 
Finally, we illustrated what in our opinion are the missing "pieces of the puzzle" required 
to understand fully the phenomenology of these remarkable materials. 
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Figure 1. Perseptive view of the crystal structure of RM^Os . Mn 3+ , Mn 4+ and R 3+ 
ions are shown respectively as red, green and pink spheres. Polyhedra around Mn 3+ 
and Mn 4+ ions connecting first neighbour oxygens ions are also shown. Magnetic 
super-exchange interactions Jj, i=3,5 are shown by double side arrows. 
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Figure 2. Neutron powder diffraction patterns as a function of temperature for 
YMn2 05 (top) and TbMi^Os (bottom) in the magnetically ordered phases. The 
scattering intensity is color coded with brighter colors representing higher intensities. 
The intensity of each diagram is normalized to the most intense nuclear reflections in 
this d-spacing range ( (0,0,1) peak at 5.8 A). The position of magnetic Bragg peaks 
are indicated by white asterix. 
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Figure 3. Dependence of the (100)-k magnetic reflection with temperature and applied 
magnetic field for TDM112O5. The scattering intensity is color coded and displayed as 
a function of d-spacing and temperature for several values of the magnetic field. The 
color scale is shown on the right. A vertical line is displayed to mark the zero-field 
ICP-CP transition temperature (Color online) 
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Figure 4. H-T phase diagram of the commensurate (CM) to incommensurate (CM) 
phase transition in TbM^Os The points have been calculated by a fitting procedure 
(see text for details). The solid line is a guide to the eyes. 
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Figure 5. Time-of- flight diffraction patterns of YM112O5 at 1.6K under magnetic 
fields between and 8 Tesla. 
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Figure 6. Representation of the commensurate magnetic structure at 25K. The model 
is derived from analysis of single-crystal neutron diffraction data. Magnetic ions are 
represented with the same color scheme defined in reference (|19i ) . Magnetic moments 
are represented by black arrows. Top panel: Magnetic configuration in the ab-planc 
(2x2 unit cells) showing the presence of nearly antifcrromagnetic zig-zag chains (see 
text for details). Several Mn 4+ -Mn 3+ -A/n 3+ -Afn 4+ fragments in the structure are 
labelled by roman numbers (la and Ha along the a-axis and lb, lib along the b-axis). 
Bottom panel: Magnetic structute projected in the ac- and be- planes. Configurations 
within each fragment (la, Ha, lb, lib) is shown for 4 unit-cells along the c-direction. 
The curled arrows at the top of each magnetic chain represent the rotation direction 
of the cyloidal modulation. For clarity, the size of the magnetic moment has been 
doubled for representing modulations in the bc-plane (lb, lib) 



REFERENCES 



31 




Figure 7. Labeling scheme and relevant interactions for the linear chain model 
discussed in this section. One of the possible spin ordering modes is shown. 




Figure 8. Magnetic site labeling scheme employed for the calculation of the DM 
energy, from (|19l )) 
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Figure 9. Phase diagram for the J\ — J2 — J 1 model for J2 = — 1 and T = (isotropic 
case). Areas in solid colors indicate commensurate phases, whereas a color gradation 
indicates incommensurability. The angle a is twice the phase offset between sites Si 
and 5*2, whereas (3 is the propagation vector q z times 27r (see text). 
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Figure 10. Left: Schematic representation of the neutron spherical polarimetry 
experiment for the two domains, here shown in the idealized case of an unpolarized 
incident beam. The real (or imaginary) component of the magnetic structure factor 
projection M fcfei rotate clockwise (counter clockwise) for Domain I (Domain II), 
creating a spin polarization of opposite signs for the scattered neutron. The direction 
of the applied electric field is also indicated. Right: Magnetic structures of YM^Os 
for the different domain configurations, projected on the ab plane (I and II) and the be 
plane (III and IV). Small arrows represent magnetic moments. The observed domain 
switching mechanism is represented by the inversion (change from the red to the blue) 
in the central chain (between configurations I and II in the ab plane). Right panel: 
Partial hysteresis loop measured on the created neutron polarization element P yx for 
the — i — j Bragg peak of an YMJI2O5 crystal, warmed to 35 K after previous cooling 
to 25 under a negative -2.2 kV/cm electric field, inset: integrated pyroelectric currents 
measured on a 0.5 mm thick YMTL2O5 crystal of the same batch on cooling to 25 K in 
a negative -2.0 kV/cm electric field (bottom/blue curve), followed by warming to 35 
K and switching to a positive +2.0 kV/cm electric field (top/red curve). The data are 
normalized to the fully saturated value at 25 K. Both hysteresis bias and the switching 
ratios are in very good agreement with the neutron data. 



